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Abstract 

Leber's Hereditary Optic Neuropathy (LHON) is one of the commonest mitochondrial diseases. It causes total blindness, and 
predominantly affects young males. For the disease to develop, it is necessary for an individual to carry one of the primary 
mtDNA mutations 11778G>A, 14484T>C or 3460OA. However these mutations are not sufficient to cause disease, and 
they do not explain the characteristic features of LHON such as the higher prevalence in males, incomplete penetrance, and 
relatively later age of onset. In order to explore the roles of nuclear encoded mitochondrial proteins in development of 
LHON, we applied a proteomic approach to samples from affected and unaffected individuals from 3 pedigrees and from 5 
unrelated controls. Two-dimensional electrophoresis followed by MS/MS analysis in the mitochondrial lysate identified 17 
proteins which were differentially expressed between LHON cases and unrelated controls, and 24 proteins which were 
differentially expressed between unaffected relatives and unrelated controls. The proteomic data were successfully 
validated by western blot analysis of 3 selected proteins. All of the proteins identified in the study were mitochondrial 
proteins and most of them were down regulated in 11778G>A mutant fibroblasts. These proteins included: subunits of 
OXPHOS enzyme complexes, proteins involved in intermediary metabolic processes, nucleoid related proteins, chaperones, 
cristae remodelling proteins and an anti-oxidant enzyme. The protein profiles of both the affected and unaffected 1 1778G> 
A carriers shared many features which differed from those of unrelated control group, revealing similar proteomic responses 
to 11778G>A mutation in both affected and unaffected individuals. Differentially expressed proteins revealed two broad 
groups: a cluster of bioenergetic pathway proteins and a cluster involved in protein quality control system. Defects in these 
systems are likely to impede the function of retinal ganglion cells, and may lead to the development of LHON in synergy 
with the primary mtDNA mutation. 



Citation: Tun AW, Chaiyarit S, Kaewsutthi S, Katanyoo W, Chuenkongkaew W, et al. (2014} Profiling the Mitochondrial Proteome of Leber's Hereditary Optic 
Neuropathy (LHON) in Thailand: Down-Regulation of Bioenergetics and Mitochondrial Protein Quality Control Pathways in Fibroblasts with the 11778G>A 
Mutation. PLoS ONE 9(9): e106779. doi:10.1371/journal.pone.0106779 

Editor: Yong-Gang Yao, Kunming Institute of Zoology, Chinese Academy of Sciences, China 
Received February 28, 2014; Accepted August 8, 2014; Published September 12, 2014 

Copyright: © 2014 Tun et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted 
use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Funding: This study was supported by The Thailand Research Fund (http://www.trf.or.th) (BRG518006 to PL and RTA5680004 to VT), Office of the Higher 
Education Commission (http://www.inter.mua.go.th/main2/index.php) and Mahidol University (http://www.mahidol.ac.th) under the National Research 
Universities initiative, and Siriraj Graduate Thesis Scholarship to AWT. The funders had no role in study design, data collection and analysis, decision to 
publish, or preparation of the manuscript. 

Competing Interests: The authors have declared that no competing interests exist. 
* Email: patcharee. ler@mahidol.ac.th (PL); thongboonkerd@dr.com (VT) 



Introduction 

Leber's hereditary optic neuropathy (LHON) [OMIM 535 000] 
is one of the commonest mitochondrial inherited diseases [1], It is 
also one of the common causes of blindness in young men, and 
more than 80% of LHON patients are male [2]. As a result of 
degeneration of retinal ganglion cell layers, the patients usually 
develop bilateral acute or sub-acute, painless loss of central vision 
[3]. 



The three missense mitochondrial DNA (mtDNA) mutations 
11778G>A, 14484T>C, and 3460G>A are the primary 
mutations responsible for 95% of LHON cases [4]. These 
mutations change amino acid R340H, M64V and A52T in 
ND4, ND6 and ND1 subunits of complex I of mitochondrial 
OXPHOS respectively. All of the LHON patients detected so far 
in Thailand carry either 1 1778G>A (>90% of cases) or 14484T> 
C [5,6] . Though presence of a primary mutation is necessary to 
develop the disease, the primary mutations per se cannot explain 
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the distinctive features of LHON [4,7]. The incomplete pene- 
trance of these mutations, and the disease's greater prevalence in 
males, imply that carriage of a primary mutation alone is not 
sufficient to cause pathogenesis, and that additional genetic and 
environmental factors are involved [4] . Level of heteroplasmy of 
mutant mtDNA, mtDNA haplogroup, nuclear DNA background 
and various environmental factors have also been reported to 
influence the clinical development of LHON [4,8-13]. Many 
previous studies have been conducted to understand these putative 
risk factors better, in particular to hunt for nuclear modifier genes. 
Nuclear genetic analyses have ranged from studies of single genes 
[14-16] to global gene expression profiling studies [17-19]. Given 
that both mitochondrial and nuclear genes encode OXPHOS 
subunits and that there is cross-talk between mitochondria and the 
nucleus, differential expression of both mitochondrial and nuclear 
genes is observed in various OXPHOS deficiency models [19]. 
The mitochondrial proteome contains approximately 1,000 
proteins, 99% of which are the products of nuclear genes [20]. 
The coordinated expressions of imported nuclear encoded 
proteins with the 13 mtDNA-encoded proteins are crucial for 
the integrity of mitochondrial function. 

In order to search for nuclear encoded mitochondrial proteins 
that may influence the development of LHON, this study 
compares the mitochondrial proteomic profiles of fibroblasts from 
affected and unaffected individuals in LHON families with those 
of unrelated controls, using 2-Dimensional Polyacrylamide Gel 
Electrophoresis (2-DE) and mass spectrometry. 

Materials and Methods 

Patients and Their Unaffected Relatives 

The samples employed in this study were cultured dermal 
fibroblasts, directiy obtained from 7 affected LHON patients from 
3 unrelated pedigrees, 3 unaffected relatives from the 3 families 
(one from each), and 5 unrelated controls. All LHON cases were 
diagnosed by an ophthalmologist (WC) and were confirmed to be 
homoplasmic carriers of the 1 1778G>A mtDNA mutation. The 3 
unaffected relatives were maternally related to the affected 
individuals and all of them were homoplasmic carriers of the 
1 1778G>A mtDNA mutation (Figure SI). As a control group, five 
individuals with no familial history of eye diseases were recruited. 
They were recruited during visits to Siriraj Hospital, Bangkok, 
Thailand, for reasons unrelated to eye diseases or chronic 
metabolic diseases. The study was approved by the Ethics 
Committee of the Mahidol University, Faculty of Medicine, 
Siriraj Hospital (No. 161/2551) and the study was conducted 
according to the principal of the World Medical Association's 
Declaration of Helsinki. Written informed consent was provided 
for all the samples used in the research. 

Fibroblast Cell Culture 

The primary dermal fibroblasts from the affected and unaffect- 
ed relatives, and the unrelated controls were cultivated and 
maintained at 37°C with 12% (v/v) fetal bovine serum (FBS) in 
Dulbecco's modified Eagle's medium supplemented with ampho- 
tericin B (1 ng/ml), penicillin (100 U/ml), streptomycin (100 ug/ 
ml), 2 mM L-Glutamine, uridine 50 Hg/ml [21] in humidified 5% 
C0 2 atmosphere at 37°C. The medium was refreshed 3 times a 
week. Cells were harvested at passage 6 for mitochondrial 
isolation. 

Immunofluorescent Staining for Fibroblast Confirmation 

For immunofluorescent staining, the cells were washed three 
times with PBS and fixed in 3.8% formaldehyde in PBS at room 



temperature for 10 minutes. After rinsing with PBS, the fixed cells 
were blocked with 1 % BSA PBS for 30 minutes. Mouse 
monoclonal anti-Fibroblast surface protein (Abeam, Cambridge, 
USA) (1:50) was incubated for 2 hours at room temperature. After 
washing with PBS for 3 times, the cells were incubated with 
secondary antibody (rabbit anti-mouse conjugated with FITC 
1:2,000 in 1% BSA PBS) and Hoechst-dye 33342 at a dilution of 
1:1,000 at room temperature in the dark. Then, the cells were 
rinsed with PBS and mounted with anti-phase solution on glass 
slide and visualized under fluorescent microscopy (Nikon 
ECLIPSE 80i, Nikon Corp.; Tokyo, Japan). 

Fibroblast Mitochondrial Isolation 

Fibroblast mitochondria were isolated by differential centrifu- 
gation [22]. Before trypsinization, the cultured cells were washed 
with chilled PBS at least four times to make sure that there was no 
residual FBS. The cultured fibroblasts were then trypsinized with 
0.25% trypsin-EDTA. Samples of 0.5 x 10 6 cells were suspended in 
1 ml of isolation buffer containing 0.25 M sucrose, 10 mM 
HEPES (pH 7.5) and 0.1 mM EDTA. Cell suspensions were 
sonicated with a probe sonicator (Bandelin Sonopuls HD 200; 
Bandelin electronic; Berlin, Germany) at MS 72/D (50 cycles) for 
10 sec. Cell lysates were then centrifuged twice at 1,000 g for 
10 min to remove cell debris and intact cells, if present. 
Supernatant was collected and centrifuged at 20,000 g for 
30 min. Pellets were saved and washed with the buffer containing 
0.25 M sucrose and 10 mM HEPES (pH 7.5) and centrifuged 
again at 20,000 g for 20 min. As a final step, the mitochondrial 
pellets were washed with PBS and centrifuged at 20,000 g for 
10 min. The whole extraction procedure was performed at 4°C. 
The mitochondrial pellets were lyzed by Laemmli or 2-DE buffer 
depending on the subsequent step of the experiment, and the lyzed 
mitochondrial fractions were kept at — 20°C until use. 

Mitochondrial Fraction Purity Evaluation 

To confirm the purity of the extracted mitochondrial fractions, 
mitochondrial proteins and the proteins from the whole cell lysate 
were resolved by western blot analysis, using mitochondrial- and 
other organelle-specific markers. The mitochondrial pellets or the 
primary fibroblasts subsequently used for the Western Blot 
experiments were lyzed by 2x Laemmli buffer containing 4% 
SDS, 10% 2-mercaptoehtanol, 20% glycerol, and 0.125 M Tris 
HC1 without bromophenol blue [23]. The protein samples were 
boiled for 5 min at 95 °C. In each western blot analysis, 20 ug of 
total proteins from each sample were loaded. The proteins were 
resolved by 3.7% SDS-PAGE (stacking gel) and 12% SDS-PAGE 
(resolving gel) at 150 V for two and a half hours by vertical gel 
electrophoresis. The resolved proteins in the gel were then electro- 
transferred to nitrocellulose membrane using a semidry transfer 
method (Bio-Rad; Hercules, CA) for 80 min with a constant 
current of 75 mA. Non-specific binding to the membrane was 
blocked with 5% skimmed milk in PBS for 1 hr. The blocked 
membranes were probed with the desired specific primary 
antibodies in 1% skimmed milk or 1% Bovine Serum Albumin 
in PBS for overnight at 4°C at concentrations according to 
manufacturers' instructions. Membranes were washed with PBS 
for 3 times (5 min per wash) and further incubated with the 
required secondary antibodies conjugated with horse radish 
peroxidase (Dako, Glostrup, Denmark), in 1 % skimmed milk in 
PBS or 1% BSA in PBS. Incubations with the secondary 
antibodies were carried out at room temperature for 1 hr, with 
the concentrations of the secondary antibodies half of those used 
for the primary antibodies. Bands were visualized by Super Signal 
West Pico chemiluminescence substrate (Pierce Biotechnology 
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Figure 1. Assessment of the purity of fibroblasts from a 
cultured skin biopsy. Fibroblast surface protein (FSP) was used as a 
marker in immunofluorescence of the cultured fibroblasts obtained 
directly obtained from the skin biopsy. The green represents fibroblasts 
and the nucleus was stained with Hoechst-dye 33342 which shows 
blue. 

doi:1 0.1 371/journal.pone.01 06779.g001 

Inc.; Rockford, IL, USA). The following primary antibodies were 
used: rabbit polyclonal anti-voltage dependent anion selective 
channel protein 1 (VDAC-1), a mitochondrial marker (Abeam, 
Cambridge, USA; ab28777); rabbit polyclonal anti-lysosomal 
associated membrane protein-2 (LAMP- 2), a lysosomal marker 
(Abeam, Cambridge, USA; ab37024); rabbit polyclonal anti- 
Calnexin, an endoplasmic reticulum marker (Abeam, Cambridge, 
USA; ab22595); mouse monoclonal anti-c-Jun, a nuclear marker 
(Santa Cruz Biotechnology, Inc, sc 166540) and mouse monoclonal 
anti-a-tubulin, a cytoplasmic marker (Santa Cruz Biotechnology, 
Inc, sc23948). 

Mitochondrial Proteomics Profile Analyses 

Two-Dimensional Electrophoresis. Mitochondrial pellets 
derived from cultured primary fibroblasts of the 7 individuals with 



LHON, the 3 unaffected relatives and the five unrelated controls 
were lyzed with a lysis buffer containing 7 M urea, 2 M thiourea, 
2% CHAPS, 120 mM DTT, 40 mM Tris, and 2% ampholyte (pH 
3—10), then incubated at 4°C for 30 min. Protein concentration 
was determined by the Bradford method [24]. Samples of 100 (a.g 
of total mitochondrial protein from each individual were mixed 
with rehydration buffer (7 M urea, 2 M thiourea, 2% CHAPS, 
120 mM DTT, 40 mM Tris-base, 2% ampholytes (pH 3-10) and 
a trace of bromophenol blue) to make the final volume of 150 ul 
The samples were rehydrated onto 7 cm immobilized pH gradient 
DryStrips (non-linear pH gradient of 3-10; GE Healthcare, 
Uppsala, Sweden) at room temperature for 1 6 hr. IPG strips were 
focused in an Ettan IPGphor II IEF System (GE health care) at 
20°C using a stepwise mode to reach 9083 Vh. 

After isoelectric focusing, the strips were equilibrated in two 
different equilibration buffers at room temperature, for 15 min 
each time. The first equilibration buffer contained 6 M urea, 
130 mM DTT, 112 mM Tris-base, 4% SDS, 30% glycerol, and 
0.002% bromophenol blue, and the second equilibration buffer 
contained the same components except 135 mM iodoacetamide 
instead of DTT. The strips were further loaded on 13% 
polyacrylamide gel and resolved using a SE260 mini Vertical 
Electrophoresis Unit (GE Health Care) at 150 V for approxi- 
mately 2 hr. The separated proteins were fixed with 10% 
methanol and 7% acetic acid for 30 min. The fixed solution was 
then removed and the gels were stained with 20 ml of Deep Purple 
fluorescence stain (GE Healthcare) overnight on a continuous 
gentle rocker. Gel images were taken using a Typhoon 9200 laser 
scanner (GE Healthcare). 

Analysis of Protein Spots. Detection and matching of spots 
on gel images and analysis of protein spots was performed using 
ImageMaster 2D Platinum software from GE Health care. 
Parameters used for spot detection were: minimal area of 10 
pixels, smooth factor of 2 and saliency of 2. The gel containing all 
of the spots and with the greatest number of spots of all the gels 
was classed as the reference gel. It was used to check for the 
presence and differential expression of proteins among gels. 
Background subtraction was performed, and the intensity volume 
of each spot was normalized with total intensity volume 
(summation of the intensity volumes obtained from all spots 
within the same 2-D gel). Intensity volumes of individual spots 




VDAC-1 (32 kDA) Mitochondria 

Calnexin (90 kDA) Endoplasmic reticulum 
C-Jun (36 kDA) Nucleus 

LAMP-2 (110 kDA) Lysosome 
a-tubulin (50 kDA) Cytoplasm 



Figure 2. Western blot analyses for assessment of mitochondrial enrichment and purity. 20 jig of mitochondrial lysate and whole cell 
lysate from fibroblasts were separated by 12% SDS-PAGE gel and checked with specific antibodies against various sub-cellular organelles. (W = whole 
cell lysate; M = mitochondrial enriched fraction) (The same membrane for each cell type was stripped and probed with subsequent antibodies.). 
doi:1 0.1 371 /journal.pone.01 06779.g002 
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Figure 3. Representative proteome map of differentially expressed proteins from the mitochondrial fractions of fibroblasts from 
(A) LHON cases (n = 7) and (C) controls (n = 5). Equal amounts of proteins (100 ug) from each fibroblast sample were resolved by 2-DE. The 
numbers indicate the spot IDs of proteins whose expression levels differ significantly between the fibroblasts of LHON cases and controls. 
doi:1 0.1 371 /journal.pone.01 06779.g003 



from each gel were subjected to statistical analysis to compare 
between the three groups of samples in the study. Protein spots 
that were differentially expressed at significant levels were 
subjected to in-gel tryptic digestion for identification by mass 
spectrometry. 

The significantly differentially expressed spots of proteins were 
manually excised from the gels. The excised gel pieces were 
washed twice with 200 uL of 50% acetonitrile (ACN)/25 mil 
NH 4 HC0 3 buffer (pH 8.0) at room temperature for 15 min, and 
then washed once with 200 U.L of 100% ACN. After washing, the 
solvent was removed and the gel pieces were dried. The dried gel 
plugs were then rehydrated with 10 |iL of 1% (w/v) trypsin in 
25 mM NH 4 HCO :i . After rehydration at 37°C for 30 min, the gel 
pieces were crushed and further incubated with 1 % (w/ v) trypsin 



at 37°C for at least 16 hr. Peptides were subsequently extracted 
twice with 50 |J.L of 50% ACN/5% trifluoroacetic acid; the 
extracted solutions were then combined and dried with the 
SpeedVac concentrator. The peptide pellets were resuspended 
with 10 |xL of 0.1% TFA and concentrated. The peptide solution 
was then washed with 10 |iL of 0.1% formic acid by drawing up 
and expelling the washing solution three times. The peptides were 
eluted with 5 uL of 75% ACN/0.1% formic acid. 

Protein Identification by Mass Spectrometry (Q.-TOF MS 
and/ or MS/MS). Digested peptides were injected into a Magic 
C18 column (Michrom Bioresources, Inc., CA, USA), which was 
connected to the MAGIC 2002 (Michrom Bioresources, Inc., CA, 
USA) high-performance liquid chromatography (HPLC) system. 
The solvent composition of the mobile phase was programmed to 
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Figure 4. Representative proteome map of differentially expressed proteins from mitochondria fraction between fibroblasts from 
(U) the unaffected LHON individuals (n = 3) and (C) the controls (n = 5). Equal amounts of proteins (100 |jg) from each sample were resolved 
by 2-DE. The numbers indicate the spot IDs of proteins whose expression levels differ significantly between the fibroblasts of unaffected LHON 
individuals and controls. Some spots in the same horizontal row showed the same proteins. For example, spot IDs 300 and 301 were identified as 
LONP1 and IDs 448, 451, 468 as catalase. 
doi:1 0.1 371 /journal.pone.01 06779.g004 
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Figure 5. Representative proteome map of differentially expressed proteins from the mitochondrial fractions of fibroblasts from 
(A) LHON cases (n = 7) and (U) unaffected LHON individuals (n = 3). Equal amounts of proteins (100 |ig) from each sample were resolved by 
2-DE. The numbers indicate the spot IDs of proteins whose expression levels differ significantly between the LHON cases and their unaffected 
relatives. 

doi:1 0.1 371 /journal.pone.01 06779.g005 



change in 50-min cycles with varying mixing ratios of solvent A 
(2% v/v CH 3 CN and 0.1% v/v HCOOH) to solvent B (90% v/v 
CH 3 CN and 0. 1 % v/v HCOOH) at flow rate of 1 uL/min using a 
MAGIC Variable Splitter. Thereafter, the peptides were eluted 
with a linear gradient from 0% to 50% solvent B. Purified peptides 
were introduced from HPLC to 0_-star (Applied Biosystems, Foster 
City, CA, USA), a hybrid quadrupole time-of-flight mass 
spectrometer, through a FortisTip (AMR, Tokyo, Japan). The 
MS/MS data were extracted and proteins were identified using 
the MASCOT search engine (http://www.matrixscience.com), 
assuming that peptides were monoisotopic, and that fixed 
modifications were carbamidomethylations of cysteine residues, 
whereas variable modifications were oxidations at methionine 
residues. Only one missed trypsin cleavage was allowed, and 
peptide mass tolerances of 50 ppm were allowed for MS/MS ion 
searches. The searches were done against human proteins in the 
NCBI database (http://www.iicbi.nlm.nih.gov). Peptides with ion 
scores greater than 34 were considered as significant hits. Only the 
significant hits from the MS /MS peptide ion search were reported. 

Western Blot Analysis for confirmation of 2-D proteomic 
results 

Western blot analysis was performed to validate the levels of 
protein expression changes identified in the proteomics profiling. 
20 ug samples of the mitochondrial fractions from each sample 
were resolved employing the same protocol described above. 
Three primary antibodies were used: rabbit polyclonal HSP60 
(Santa Cruz Biotechnology, Inc, sc 13966), rabbit polyclonal anti- 
Catalase (Abeam, Cambridge, USA; abl6731) and rabbit poly- 
clonal NDUFS1 (Abeam, Cambridge, USA; ab96428). Rabbit 
polyclonal anti-VDAC was used as a loading control. After 
incubating with the respective secondary antibodies (anti-rabbit), 
bands were visualized by enhanced chemiluminescence and 
exposed to film. The band intensities were measured using ImageJ 
software (http:/ /rsbweb.nih.gov/ij/). 



Statistical Analysis 

All data representing the intensity volume of the spots were 
reported as mean ±SEM. For comparison between 3 different 
groups, the data were analyzed using one-way ANOVA followed 
by a Post Hoc Tukey- Kramer Test (SPSS, version 18). P- values 
less than 0.05 were considered statistically significant. The false 
discovery rate (FDR) was determined using R statistical packages 
for calculating q-values and FDR-values [25,26]. 

Results 

LHON Patients and Their Unaffected Relatives 

Samples from individuals in three unrelated LHON pedigrees 
with different mtDNA haplogroup backgrounds were employed in 
this study (Figure SI). The numbers of affected LHON patients 
recruited from pedigrees Fl, F9 and F66 were two, two and three 
respectively. All affected individuals have been diagnosed by an 
ophthalmologist (WC) and have been confirmed as bearing 
homoplasmic mtDNA 117780A [9,13,27]. The characteristics 
of each patient are summarized in Table 1. The 3 unaffected 
relatives came from the 3 pedigrees and were maternally related to 
the affected individuals and also confirmed to be homoplasmic 
carriers of the mtDNA 117780A mutation. The 11778G>A 
mutation was not detected in the 5 control individuals recruited in 
this study, after performing both RFLP genotyping and Sanger 
sequencing in both the forward and reverse directions. 

Confirmation of the purity of Fibroblasts 

One of the fibroblast cultures from each of the skin biopsies was 
randomly selected and tested with anti-fibroblast surface protein 
which is specific to fibroblasts [28]. Figure 1 shows the results of 
this immunofluorescent staining with DNA staining Hoechst-dye 
33342, which indicated that almost every stained cell had a 
positive signal for fibroblast surface protein. This confirmed that 
the fibroblast cultures were pure and uncontaminated by other cell 
types. 
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Mitochondrial enrichment and purity 

Enrichment and purity of the mitochondrial fraction was 
assessed by immunoblotting, using markers specific for mitochon- 
dria and other organelles. Figure 2 demonstrates that the 
mitochondrial marker VDAC-1 was highly enriched in the 
mitochondrial fraction compared with the whole primary fibro- 
blast lysates. In addition, the markers for endoplasmic reticulum, 
nucleus, lysosome and cytosol were absent or minimal in the 
mitochondrial enriched fraction. The mitochondrial fraction was 
highly enriched and with minimal level of non-mitochondrial 
contamination, and hence it was suitable for further analysis with 
2-DE. 

Mitochondrial Proteomics Profile Analyses 

2-D PAGE comparisons of mitochondrial fraction 
between fibroblasts from LHON cases, unaffected 
relatives and unrelated controls. The proteins from the 
mitochondrial enriched fraction of the primary fibroblasts were 
analyzed by 2-DE using pH 3-1 1 non-linear pH gradient strips for 
the first dimension and 13% SDS-PAGE for the second 
dimension. The individual 2-D mitochondrial proteome profiles 
from fibroblasts of individuals with LHON (n = 7), unaffected 
relatives (n = 3) and unrelated controls (n = 5) showed virtually 
identical spot patterns (Figure 3, 4 and 5). Approximately 800 
protein spots were visualized on each gel. However, spot matching 
analysis revealed distinct differences in the proteomic profiles of 
the three groups of samples. 

Interestingly, the majority of all differential-intensity spots were 
lower intensity in the LHON cases than in the unrelated controls. 
The degrees of fold-change between the cases and the unrelated 
controls ranged from 0.47 to 2.05. Meanwhile, the fold-changes 
between the unaffected relatives and the unrelated controls ranged 
from 0.28 to 3.90, and those between cases and unaffected 
relatives ranged from 1.5 to 2.2 fold. In total 61 differential 
intensity spots were identified in the three pairwise comparisons. 
These were excised from the gels, and their corresponding proteins 
were identified by MS or MS/MS analysis. 

Proteins which were differentially expressed in the 
fibroblast of LHON cases, unaffected relatives and 
unrelated controls. Some sets of spots identified by MS/MS 
analysis in the same horizontal row were found to correspond to 
the same protein. Of the 61 differential-intensity spots, 33 spots 
representing 27 unique proteins were identified. Different isoforms 
of the protein mitofilin were represented by different spots in the 
gel. Seventeen proteins were differentially expressed between 
LHON cases and unrelated controls, and 24 proteins were 
differentially expressed between unaffected relatives of cases and 
unrelated controls (Table 2 and 3). Seven proteins were differently 
expressed between the LHON cases and their unaffected relatives 
(Table 4). An assessment of mitochondrial localization using the 
MitoMiner database [29] indicated that most of these identified 
proteins were mitochondrial resident proteins, and some had 
functional associations with mitochondria, though some were also 
localized in other compartments of the cell (Table 5). 

Functions of the proteins were assessed using Nextprot [30]. 
The differentially expressed proteins could be classified in several 
groups: (1) those involved in intermediary metabolism (of 
carbohydrates, lipids and amino acids); (2) subunits of OXPHOS; 
(3) a cristae remodeling protein; (4) a protein involved in 
mitochondrial gene expression; (5) signal transduction; (6) chaper- 
onins and proteases of the protein quality control system; (7) an 
antioxidant enzyme and (8) cytoskeletal structure (Table 5). All of 
these proteins were down-regulated in the mitochondrial fibroblast 



11 



September 2014 | Volume 9 | Issue 9 | e1 06779 



Mitochondrial Proteomics in LHON 



Q 



a 
u 
at 
8: 
< 

°>TJ 

.c t5 
u ,2 

2"S 
o c 

Li. 3 



3 o 



0) 
T3 

tl 



T3 



O 
X 



0) 

o 

Q. 

T3 

OJ 

CD 

To 
_>~ 

c 



C 



O 

>s 

rc 
E 

E 
c/i 

-Q 

IS 



5 ■ -S S 

£ *! a i/i 

3 = ■ S 

c £ is 



> s 

u g 



a 5 



U 
z 



o 










LTl 




o 


o 




OJ 


un 


O 


m 


O 


O 


o 


O 


O 


o 


d 


+ 1 


+1 


+ 1 


m 


rsi 


a\ 




r-» 




rs 


LO 




o 






o 


d 


d 



Q 
Q 



< 
X 

a 



"D 
I 



< 
> 

o 



cu 



^ 03 — 



(T3 'o> 



O g 

'ro 1 O 
£ a. 



y a) -§ cU y - 
.9- en >, a> □_ .E 

O 2 £ 2 oi r 

■n "5 Q "9 c | 

5 



O 5 Si 

X D C 



^ x: -a .c 
r id c cj 

a a < o 



CO 



cu t a; 

£ c s= ® 

o tO p Q. 

£ E £ E 

3 d a» o 

in I Q U 



rc O 

£ -a £ 

n: >, 

H -C o 

2 ci a 

<J D < 



C 4-1 

si 



- 



a o & : 



« ■= ™ -n £ 



3 C? £ a 3 



o -3 



£ 5 S 

13 1 " 

0 »N 

§ £ £ 

1 £? 
° ■£ ^ 

4-1 U (U 

C ra S 

3 -° 3 

0 S c 

y « =3 

ro > o 

£ t/i 

« Jn 

01 U o 

£ 1 s 



PLOS ONE | www.plosone.org 



12 



September 2014 | Volume 9 | Issue 9 | e1 06779 



Mitochondrial Proteomics in LHON 



Table 5. Functional categories and sub-cellular localizations of proteins 
and Nextprot. 


identified in the study based on the MitoMiner Database 




Proteins 


Localization 


Intermediary metabolism: TCA cycle and Carbohydrate Metabolism 


2-oxoglutarate dehydrogenase, mitochondrial isoform 1 precursor 


M 


Glycerol-3-phosphate dehydrogenase 


M 


Dihydrolipoyl dehydrogenase 


M 


Pyruvate dehydrogenase El component subunit alpha 


M 


Intermediary metabolism: Fatty acid Catabolism 


Methylmalonyl-CoA mutase 


M 


Trifunctional enzyme subunit alpha 


M 


Very long-chain specific acyl-CoA dehydrogenase 


M 


Propionyl-CoA carboxylase 


M 


Intermediary metabolism: Amino acid Metabolism 


Glutaminase 


M 


Glutamate dehydrogenase 1 


M 


Subunits of oxidative phosphorylation and electron transport function 


Succinate dehydrogenase [ubiquinone] flavoprotein subunit 


M 


Cytochrome b-cl complex subunit 1 


M 


Electron transfer flavoprotein subunit alpha 


M 


ATP synthase subunit alpha 


M 


NADH dehydrogenase (ubiquinone) Fe-S protein 1, 75 kDa {NADH-coenzyme Q reductase} 


M 


Cristae remodeling 


Mitochondrial inner membrane protein 


M 


Mitochondrial Gene expression 


Leucine-rich PPR motif-containing protein 


M 


signal transduction 


Major vault protein 


N, M 


Glucosidase 2 subunit beta 


ER 


Protein stability and degradation of protein 


Lon protease homolog 


M 


Stress-70 protein, mitochondrial/Heat shock 70 kDa protein 9/MTHSP75 


M, C 


60 kDa heat shock protein, mitochondrial 


M, C 


Anti-oxidant enzymes 


Catalase 


P, Ml 


cytoskeletal protein 


Vinculin 


CSK 


(5-actin 


CSK 


Vimentin 


CSK 


Others 


Myosin-9/Cellular myosin heafvy chain-type A 


C, M 


M = Mitochondria, Ml = Mitochondrial intermembrane space, ER = Endoplasmic Reticulum, C 
doi:1 0.1 371 /joumal.pone.01 06779.T.005 


= Cytoplasm, N = Nucleus, CSK = cytoskeleton. 



proteomes of LHON patients and their unaffected relatives, 
compared with the unrelated controls. 

However, expression levels of heat-shock protein 60, methyl- 
malonyl-CoA mutase and myosin heavy chain were significandy 
different between the cases and the unrelated controls, but not 
between the unaffected relatives and the unrelated controls. 
Conversely, expression levels of 2-oxoglutarate dehydrogenase, 
NADH dehydrogenase (ubiquinone) Fe-S protein 1 (NDUFS1), 
ATP synthase subunit alpha, glutaminase, glutamate dehydroge- 
nase, propionyl CoA carboxylase, beta-actin, glucosidase 2 subunit 



beta (PRKCSH) and vimentin were significantly different only 
between the unaffected relatives and the unrelated controls. 

Expression of seven proteins differed significantly between the 
cases and their unaffected relatives: these were subunits of the 
respiratory chain such as ubiquinol cytochrome c reductase core I 
protein and ATP synthase subunits alpha; an enzyme of the TCA 
cycle, dihydrolipoamide dehydrogenase; one of the subunits of the 
pyruvate dehydrogenase complex; lon protease and major vault 
protein (Table 4). 
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Control LHON 




60 (60 kDa) 



Figure 6. Validation of proteomic data by Western Blot analysis. These figures show: (A) The differentially expression of heat-shock protein 
60 and catalase between LHON cases and controls; and (B) the differential expression of NDUFS1 and catalase between unaffected mutation carriers 
and controls. VDAC was used as a loading control. The two lanes represent samples from two different individuals for each group of samples. 
doi:1 0.1 371 /journal.pone.01 06779.g006 



Validation of Proteomic data by Western Blot 

To validate the differentially expressed proteins observed in the 
2-DE proteomic analysis, three of the differentially expressed 
proteins - heat shock protein 60, catalase and NDUFS1 - were 
selected for western blot analysis. Compared with the unrelated 
controls, heat shock protein 60 and catalase were down-regulated 
in the mitochondrial fraction of the LHON cases, while catalase 
and NDUFS1 were down-regulated in the unaffected relatives. 
These differences all ran in the same direction as in the 2-DE 
proteomic analysis (Figure 6). 

Discussion 

It is impossible to sample retinal ganglion cells, the affected 
target cells of this mitochondrial disease. Hence in this study we 
used skin fibroblasts derived from the patients, which may mirror 
some of the features of other post-mitotic tissues such as retinal 
ganglial cells. The fibroblasts were cultured directly from the skin 
biopsies, from LHON patients and their unaffected relatives 
homoplasmic for the primary mtDNA mutation 1 1 778G>A. The 
purity of the skin fibroblast cultures, which may contain 
epithelium, reticulocytes, and adipocytes, was confirmed using 
anti-Fibroblast surface protein. The skin fibroblasts were then used 
to explore the mitochondrial proteomes of these 11778G>A 
carriers, comparing them with the proteomes of control individ- 
uals using 2-DE with a non-linear IPG strip of pH 3-1 1. 

The purity and the enrichment of the mitochondrial fraction 
were confirmed by western blotting with organelle-specific 
markers. The proteomics data also confirmed this, as most of 
the identified proteins were mitochondrial resident proteins. No 
mtDNA-encoded proteins, including PARL (hydrophobic and 
membrane bound protein) which was previously reported [27] 
were identified in this study, since these proteins are highly 
hydrophobic [31-33] and 2-DE does not provide good separation 
of hydrophobic and membrane bound proteins [34-36] . 



All of the proteins that were differentially expressed between 
11778G>A carriers and unrelated controls were down-regulated 
in fibroblasts of 11778G>A carriers: either the LHON cases or 
their unaffected relatives (Figure 7). Some of the subunits of the 
OXPHOS complex were altered: NDUFS1 in complex I, 
succinate dehydrogenase in complex II, UQCRC1 in complex 
III and ATP5A1 in complex V were all down-regulated in 
fibroblasts of either the affected or the unaffected 11778G>A 
carriers. This was an unexpected finding since complex I defective 
LHON mutant cells may compensate for a lack of ATP synthesis 
through up-regulation of succinate/glycerol-3-phosphate dehydro- 
genase [37-39]. However, our results were in agreement with the 
findings of Qiang et al, where succinate/glycerol-3-phosphate 
driven respiration was reduced in mutant cells derived from some 
Chinese LHON families [40]. 

The proteins found to be significandy differently expressed 
between the cases, unaffected 11778G>A carriers and unrelated 
controls can be categorized functionally into 2 functional groups: 
bioenergetic pathways and mitochondrial protein quality control. 
A number of enzymes for intermediary metabolism were 
differentially expressed in 11778G>A mutant fibroblasts com- 
pared with unrelated controls (Table 2 and 3). Interestingly, 
differential expression was observed particularly in FAD-linked 
and NAD-linked dehydrogenases which produce reducing equiv- 
alents to feed the respiratory chain. For instance, altered 
expression was found for two enzymes of fS-oxidation: FAD- 
dependent very long chain specific acyl-CoA dehydrogenase, and 
NAD-dependent trifunctional enzyme subunit alpha (also known 
as hydroxyacyl CoA dehydrogenase). And among the TCA cycle 
enzymes, FAD-dependent succinate dehydrogenase, NAD-depen- 
dent 2-oxoglutarate dehydrogenase, pyruvate dehydrogenase El 
component alpha subunit and dihydrolipoyl dehydrogenase were 
down-regulated in 11778G>A fibroblasts. Also down-regulated 
was glycerol-3-phosphate dehydrogenase, which is involved in 
transporting cytosolic reducing equivalents to the mitochondrial 
respiratory chain. 
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Figure 7. (A) A plot of log-ratio spot intensities, comparing affected individuals (black bars) and unaffected relatives (grey bars) 
versus controls. Log-ratios of protein spot intensities were plotted according to the method of [63]. In the case that one protein was identified in 
several spots, the only spot intensity that gave the highest significant value was selected for this plot. Negative changes were observed in all 
discoveries. Group A includes protein spots which were significantly differentially expressed in comparisons of affected and unaffected carriers, 
affected individuals and controls, and unaffected individuals and controls (Post Hoc Tukey Test; P-value<0.05). Group B includes protein spots which 
were significantly differentially expressed in comparisons of affected and unaffected carriers, and between unaffected carriers and controls. Changes 
between affected individuals and controls were not significant. Group C includes protein spots which were significantly differentially expressed only 
in the comparison between affected individuals and control. Group D includes protein spots which were significantly differentially expressed only in 
the comparison between unaffected carriers and controls. Group E includes protein spots which were significantly differentially expressed the 
comparisons between affected individuals and controls, and between unaffected carriers and controls. Here, the differences in expression between 
affected and unaffected carriers are not significantly different. (B) Venn diagram representation of proteins identified in three comparisons 
derived from figure 7A. 
doi:1 0.1 371 /journal.pone.01 06779.g007 



Given the decreased levels of metabolic and OXPHOS related 
proteins in mutation carriers in the present study, it is likely that 
aerobic respiration is particularly affected in LHON mutant 
fibroblasts. To further support this notion, ETFA electron transfer 
flavoprotein subunit alpha was down regulated in the fibroblasts of 
both affected and unaffected mutation carriers. This protein is 
important in transferring electrons from many of the mitochon- 
drial flavoprotein dehydrogenases to the respiratory chain [41]. 
This bioenergetic dysfunction is consistent with previous reports of 
defective energy metabolism in LHON mutant fibroblasts, 
including reduction of complex I activity, poor respiratory 
capacity and reduced ATP content in mutants compared to 
controls [42,43]. 



Several stability and transport proteins were also differentially 
expressed: heat shock protein 60 (HSPD1), Stress-70 protein 
(MTHSP75/HSPA9) and Ion protease 1 (LONP1) (Table 2 and 
3). HSPD1 is one of the most important chaperonins inside the 
mitochondrial matrix. It facilitates correct folding and prevents 
mis-folding of unfolded proteins formed under mitochondrial 
stress [44]. Mutations of HSPD1 in hereditary spastic paraplegia 
and MitChap60 disease highlight possible implications for 
neurodegenerative disease [45-47]. In the present study, expres- 
sion of HSPD1 was reduced 1.5 fold in LHON case fibroblasts 
compared to unrelated controls. This reduction may have 
deleterious consequences if it happens in neuronal retinal ganglion 
cells, since they are post-mitotic and are highly susceptible to the 
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Table 6. Differentially expressed proteins which show similar trends of down-regulation in the mitochondrial proteomes of both 
LHON cases and unaffected relatives when compared with controls. 



Name of Protein Cellular pathway or function involved 



Glycerol-3-phosphate Intermediary metabolism: transferring reducing equivalents to mitochondrial repiratory chain 

dehydrogenase 



Dihydrolipoyl dehydrogenase 


Intermediary metabolism: TCA cycle; generation of reducing equivalents for OXPHOS 


Pyruvate dehydrogenase 
E1 component subunit alpha 


Intermediary metabolism: carbohydrate; generation of reducing equivalents for OXPHOS 


Trifunctional enzyme 
subunit alpha 


Intermediary metabolism: lipid; p-oxidation; generation of reducing equivalents for OXPHOS 


Very long-chain specific 
acyl-CoA dehydrogenase 


Intermediary metabolism: lipid; p-oxidation; generation of reducing equivalents for OXPHOS 


Succinate dehydrogenase 
[ubiquinone] flavoprotein subunit 


Intermediary metabolism: TCA cycle; generation of reducing equivalents for OXPHOS 


Electron transfer flavoprotein 
subunit alpha 


Transfer electron from flavoprotein to OXPHOS 


Mitofilin 


Cristae remodeling 


Lon protease 1 


Protein stability and degradation 


Stress-70 protein/Heat shock 
70 kDa 

protein 9/MTHSP75 


Protein stability and degradation 


Catalase 


Anti-oxidant enzyme 


Leucine-rich PPR motif- 
containing protein 


Mitochondrial gene expression 


Vinculin 


Cytoskelon 


Major vault protein 


Signal transduction 
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accumulation of unfolded proteins [48]. Another chaperone 
Stress-70 protein (MTHSP75/HSPA9) was down-regulated in 
fibroblasts of both affected and unaffected mutation carriers 
compared with controls. In addition, HSPA1, HSPA9 and 
LONP1 co-occurred as differentially expressed proteins in other 
cellular models [49] . They are the components of the mitochon- 
drial protein quality control system which protects the formation 
of the unfolded protein with chaperones and clears protein 
aggregates with proteases [50]. There can be detrimental 
consequences when the system is under-expressed, as found in 
the 11778G>A fibroblasts of the present study. 

In addition to the proteomic changes in metabolic enzymes and 
the protein quality control system, some proteins controlling 
mitochondrial gene expression were also down-regulated in the 
mutation carriers: Leucine-rich pentatricopeptide repeat motif- 
containing protein (LRPPRC) and LONP1. LRPPRG was down- 
regulated in the fibroblasts of both affected and unaffected carriers 
of 1 1778G>A. It is a disease modifier in Leigh syndrome French- 
Canadian type [51]. Silencing of LRPPRC was associated with 
reduction of mitochondrial proteins including mitochondrial and 
nuclear encoded subunits of OXPHOS [52]. LONP1, a nucleoid 
protein, was also under-expressed in mutant cells. LONP1 binds to 
mtDNA, and its level influences a cell's sensitivity to mtDNA 
damage [53], with potential implications for LHON pathogenesis. 

Apart from proteins involved in bioenergetic pathways and 
mitochondrial protein quality control, catalase and mitofilin were 
also differentially expressed. Compared to controls, catalase was 
down-regulated in both affected and unaffected fibroblasts 
carrying 117780A. Though it is mainly located in the 
peroxisome [54], it is also associated with mitochondria in some 
cell types [55,56]. Reduction or absence of catalase under 
oxidative stress in mitochondria can lead to the inefficient 



degradation of H 2 0 2 , leading to potential mitochondrial and 
cellular damage [57]. 

Mitofilin, a protein which plays a role in/influences cristae 
morphology and nucleoid structure, was down-regulated. Previ- 
ously, down-regulation of mitofilin was observed in dopamine 
induced oxidative stress [58] and MPTP induced complex I 
inhibition [59] . Since oxidative stress and complex I inhibition are 
associated with LHON mutations, the down-regulation of mitofilin 
in 117780A mutant cells may have a role in LHON 
pathogenesis. 

It is interesting to observe the similarities in fibroblast 
mitochondrial protein expression profiles between affected and 
unaffected 1 1778G>A carriers, in comparison with the unrelated 
controls. The levels of 14 (out of total 27) identified proteins were 
reduced both in the affected and the unaffected mutation carriers 
compared to the controls (Table 6). The proteins down-regulated 
in both groups of 1 1778 mutant fibroblasts were electron transfer 
flavoprotein, succinate dehydrogenase, dihydrolipoyl dehydroge- 
nase, subunit of pyruvate dehydrogenase, glycerol-3-phosphate 
dehydrogenase concerned with electron transfer and aerobic 
respiration, very long chain specific acyl-CoA dehydrogenase, 
trifunctional plasma enzyme of P-oxidation, chaperonins such as 
lon protease 1, heat-shock protein-70, leucine-rich PPR motif- 
containing protein, mitofilin, catalase and some other proteins 
(Table 6). This similarity might reflect common compensatory 
responses cells responding similarly to the same adverse condi- 
tions, regardless of whether the cells come from affected or 
unaffected individuals. Since 11778G>A mutant fibroblasts have 
been shown to have reduced complex I activity and larger 
bioenergetic defects [43], these reduced levels of protein expres- 
sion might be due to degradation of proteins with poor 
performance or to alteration of nuclear gene responses in the 
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presence of OXPHOS deficiency. These mechanisms would be 
consistent with the observed results: that the cells in primary 
fibroblasts cultures encountering the 11778G>A mutation 
respond similarly to the adverse consequences of the mutation, 
regardless of the phenotypic status of the person who provided the 
biopsy. This mutation might reprogram gene expression profiles, 
with similar consequent physiological states of the cells driving 
similar proteome changes between the carriers. Similar transcrip- 
tomic profiles have been observed in previous studies, where 
clinically affected and unaffected tissues from the same individual 
have similar profiles [60] . 

There are some previous studies on transcriptomal changes in 
LHON mutant cells using various cell types. However, they 
reported differing patterns of transcriptomal changes, possibly 
reflecting differences in susceptibility to bioenergetic derangements 
between cell types. For instance, transcriptomic changes from 
osteosarcoma-derived LHON cybrids and LHON-mutant lym- 
phoblastoid cell lines generally showed different profiles, with the 
exception of nine common genes [17,19]. Differentially expressed 
gene products involved in transcription and transport processes 
were reported in the lymphocytes of four Saudi Arab LHON 
patients [18]. None of these gene products was differentially 
expressed in our data. This could be due the employment of 
different cell types (fibroblasts versus osteosarcoma derived cybrids 
or lymphoblastoid cell lines), different sub-cellular locations that 
we are observing (enrichment of mitochondrial proteins in our 
study versus global expression profiling), or different patterns of 
expression at the RNA level and the protein level, because of post- 
translational modifications, degradation and dependence on 
organellar transport. Consequently, we did not find any evidence 
of ER stress or protein unfolding responses, which are observed in 
various mitochondrial diseases [19]. 

This study which, to the best of our knowledge, is the first study 
to profile mitochondrial proteomes in LHON, has some limita- 
tions. One potential limitation is the challenge of selecting well 
age-and sex-matched controls for cases and their related unaffect- 
ed individuals. We selected pedigrees in which the unaffected 
individuals were well above the age of onset of their affected 
relatives for each particular family (Table 1). All of the three 
unaffected relatives, U 1 from pedigree F 1 (30 years old at present), 
U2 from pedigree F9 (2 1 years old at present) and U3 (49 years old 
at present) from pedigree F66 were above both the mean age of 
onset for Thai individuals (20. 7± 10.0 years for males and 
28. 6± 14.6 years for females) [61], and the mean age of onset 
within their families, 13 years, 17 years and 32 years respectively. 
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